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at LHC in the near future. 
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1 Introduction 



The heavy pseudoscalar meson B c contains two heavy quarks of different 
flavor. This meson has been discovered in 1998 via the decay mode B c — ► 
J/ipl^u in 1.8 TeV pp collisions, using the CDF detector at the Fermi Lab 
PQ. The B c meson constitutes a very rich laboratory for studying various 
decay channels. There are three classes of decays of B c meson, namely b- 
quark decay (when c is spectator), c-quark decay (when b is spectator) and 
the weak annihilation channels. Because of two heavy quark contents, the 
B c decay channels are expected to be very rich in comparison with other 
B mesons, so investigation of this meson is essential from both theoretical 
and experimental point of view. The B c meson decays provide windows for 
reliable determination of the CKM matrix element V c b and can shed light on 
new physics beyond the standard model. 

At LHC with the luminosity values of C = 10 cm -2 s -1 and y/s = 14TeV, 
the number of Bf events is expected to be about 10 8 ~ 10 10 per year, so 
there are high probability to study not only some rare B c decays, but also CP 
violation, T violation and polarization asymmetries. Some possible channels 
are B c -> luj, B c -> p+ 7 , B c -> K* + j and B c -> B*l + l~ } B c -> B* ul , which 
have been studied in the frame of light-cone QCD and three point QCD sum 
rules [21 El HJ [13] . A large set of exclusive nonleptonic and semileptonic decays 
of the B c meson, which have been studied within a relativistic constituent 
quark model can be found in [5]. Another possible decay channel of B c 
is B~ — ► r(i~v decay, which is studied both for decay rate and lepton 
polarization asymmetry [6]. We analyzed the radiative B c — > D*^ decay by 
using QCD sum rules method. Using our calculations at the end, we also 
analyze the B c — > D*j decay by making the necessary changes. 

The main quantities in analyzing of B c — > D*^ decay are the form factors. 
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For the calculation of form factors, relevant to this transition, we need some 
nonperturbative approaches. Among the nonperturbative approaches, QCD 
sum rules method has received special attention, because this approach is 
based on QCD lagrangian. This method has been successfully applied to a 
wide variety of problems (for a review see [2l [3l SI d El [91 [10l HH H2] ) . 

The B c — > .0*7 decay occurs via flavor changing neutral current (FCNC) 
transition (b — > 57) and weak annihilation channels. The b quark decay 
(electromagnetic penguin) for B c — > B* / y has been calculated in [13] (for 
more details about the electromagnetic penguin diagram see also [II]). We 
repeated the similar calculations for our problem and found that the corre- 
sponding branching ratio contribution was less important than that of the 
weak annihilation channel. Note that, the B c — > 0*7 decay has been in- 
vestigated in the perturbative QCD (PQCD) approach [T5] , relativistic in- 
dependent quark model (RIQM) [16] and in the framework of pertubative 
QCD in SM (PQCD), multiscale walking technicolor (MWTCM) and top- 
color assisted MWTCM (TAMWTCM) models [Hj. They found also that 
the contribution of the weak annihilation was important than that of the 
electromagnetic penguin in PQCD, RIQM and TAMWTCM. In addition, 
B c — > 0*7 decay has also been investigated in the relativistic independent 
quark model (RIQM) [16]. 

The paper is organized as follows: In section (2) we construct the tran- 
sition amplitude for the weak annihilation channel in terms of four relevant 
form factors, where a photon can be radiated from B c or D*. Two of the 
relevant form factors for this decay are calculated in [2] in the framework of 
the light-cone QCD sum rules. In section (3), we calculate the remaining two 
form factors, when a photon is radiated from D* meson also in light-cone 
QCD. In section (4), we calculate the transition form factors for electro- 
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Figure 1: The weak annihilation mechanism for B c — > D*^ 

magnetic penguin in the framework of three point QCD sum rules method. 
Finally, in section (5) numerical analysis, discussion and comparison of our 
results to those of the other approaches are forwarded and conclusion is pre- 
sented. 

2 Transition amplitude of the weak annihila- 
tion for the radiative B c — > D*~f decay 

In this section, we concentrated on the main points for obtaining the matrix 
elements of the radiative decay B c — > D*^ along the lines similar to [18J. 
The weak annihilation mechanism for this decay is shown in (Fig.l). The 
transition amplitude for this decay can be written as: 

M(B C -> D* sl ) = < D* s (p)j(q) | (sT u c)(cT v b) \ B c (p + q) > (1) 

where T u = j v (l — 75) and p, q and p+q are the momentum of D*, photon 
and B c , respectively. Using factorization hypothesis, The matrix element in 
Eq. (pQ) can be written in the following form: 

<D:(p) 1 (q)\(sT u c)(cT"b)\B c ( P + q)> = - e e^ D ^ f D; m D *T^ 
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-iee»(p + qyf Bc TW (2 ) 

where, the covariant decomposition of hadronic matrix elements Tf£ c > and 
TjP*> are responsible for the emission of photon from initial and final states, 
Id* i fn c are the leptonic decay constants of D* and B c mesons, respectively, 
and and e^P*)" are the polarization vectors of a photon and D* meson. 
The covariant decomposition of hadronic matrix elements are defined by the 
following two-point correlation functions: 

T^ip, q)=ij d'xe^ < | T {j^cT u b(0)} I B c (p + q) > (3) 

and 

T™(p, q) = if d'xe^ < D* s {p) \ T {f™sT „c(0)} | > (4) 

where jff 1 stands for electromagnetic current. Our aim is to construct the 
T^ c ^ and Tj^^ in terms of form factors and other physical quantities. Let 
first focus on T^ c \ This quantity can be written in terms of two independent 
4- moment a p and q in general as follows: 

T jxv c) (P, q) = 9nv a + Vy.qJ> + qtiPuC + p^p v d + q^q v e + e Uf , Xa p x q a F^ Bc) (5) 

where a, b, c, d, e and Fy 3 ^ are invariant amplitudes. Applying the Ward 
identity for electromagnetic current to Eq. (J5J) and using the fact that for a 
real photon q 2 = 0, we rewrite Eq.(!51) in the following form: 

T iu c) (P, q) = {g^ip-q) - p^iF^ + g^(p.q)a + p^q v f3 + q^q u c 

+i^f Bc + q»q»e + e v ^p\ a F^ (6) 
p.q 

where _F\ , a and (3 are the new invariant amplitudes. To obtain the relation 
between a and /?, we compare Eqs. (jSJ) and (jSD, which leads to 

9^a + p^qj) = (g^(p.q) - Py,q v )iF^ c) + g^{p.q)at + p^fi (7) 
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by multiplying both sides of Eq. (0) with q^ and using a + (p.q)b = ifs c we 
get the following relation between a and (3, which is called the Ward identity. 



From Eq. (jSJ), it is clear that a and /3 can take different choices. Within the 

'fee 

(p-q) ' 



scop of the present work, in parallel with [18], we set (3 = and a = -f^k- 



Substituting the values of a and /3, we obtain 

T !iu c) (p, q) = (g^uip-q) - VvA^a + l fB c g^ + q^ c 

+i^f Bc + q,q v e + e vlxX(7 p x q° F^ (9) 
p.q 

Using e.q = 0, e^ D *\p = and Eq. we obtain the following expression 
for the first term in Eq. (j2j) in terms of two form factors (Fj^ c \ Fy^) : 

^ f Dt m D *T^ = ef Dt m Dt {[{e.e^){p.q) - (e.p)(e^ .q)]iF^ 

+if Bc (e.eW) + e^V^'} (10) 

Omitting the details for the calculation TjP°\ we get the following result for 
the second term in the Eq. (j2J): 



iee»{p + qyf B T^ = ief B M^ (Dt) )(p-q)-(e.p)(e^.q)]iF^ 



(D* s ) 



+ fD*m Dt {e.e^) + e^e^V/^ 1 } (11) 

where, and Fy are two form factors of £)*. Now, we can write the 

transition amplitude for the radiative B c — > D*7 decay in terms of four form 
factors c) , F^ c) , F { f t] and as follows: 

M(B C -> d; 7 ) = e ^^[-/ DsmD .{[( e . e W))(p. g ) - (£.p)(e^).g)]iFf c) 

+i/B e (e.e^) + e^^V^ 1 } 
-i/B c {[(£.£ ( ° ;) )(p.g) - (e.p)(e (D!) .g)]iFf :) 
+/ D: m J ,.(e.e^) + e^e^e^q^F^}} (12) 



The form factors F^* and F V B ^ corresponding to the emission of the 
photon from b and c quark (see Fig. 1 i, ii), are calculated in [2J. Therefore, 
we will concentrate on the calculation of the form factors F^ 3 ^ and Fy^"^ 
(see Fig.l iii, iv). 



3 Light cone QCD sum rules for the form fac- 
tors F { A Dt) and F^ s) 

Based on the general idea on QCD sum rules method, we will calculate the 
transition form factors by equating the representation of a suitable correla- 
tor in hadronic and quark-gluon languages. For this aim, we consider the 
following correlation function: 

lWp,g) = i J d*xe lQx < 7 (g) | T{c(x)^(l - 7 5 )s(x)s(0) 7 ,c(0)} | > 

(13) 

Here, q and p are the momentum values of photon and D*, respectively and 
Q = p + q is the transferred momentum. Now, we insert the hadronic state 
D*(p) to Eq. (|T3|) . This can be re- written as: 

n^, q) = < M I gr»(i - *>)> I p:(p) >< £M I ^> 1 > . (M) 

m D , — p l 

The second term in Eq. (|14|) . by definition is: 

< D* s (p) | s lv c | >= f Dt m Dt ef^ (15) 

From Fig. 1 (iii, iv) and due to the fact that parity, Lorentz and gauge 
invariance are musts. We can write the matrix element for the emission of 
the photon from D* meson as: 

< 1 (q)\c^(l- l5 )s\D* s (p)> = efe^e^V ^ 



IK 



(Q 2 



m 



Substituting Eqs. (TT5|) and (TToT) to (TBI) , we have 

rr i \ e fD*rriD* f . a m*\ r D *\ a Fy >s (Q 2 



mjj* — p* " m|,, 



+ [e M 4 D:) (e (2?:) -9) - (e.e ra )4 D, \]^T^} (17) 



i 

mi,. 

Summation over polarization of -D* meson is performed by using: 

«W = -W + £P (18) 

^ s 

After performing the standard calculations for the phenomeno logical part, 
we get: 

n^(p,g) = — 5-* ^{le^atrS Q 2 m E » ~ 2 J ( 19 ) 

m™ — p z m n » nip,, 

The theoretical part (QCD side) of the correlator is calculated by means 
of OPE up to operators having dimension d = 5 in deep Euclidean space, 
where both p 2 and Q 2 are large and negative. It is determined by the bare- 
loop (fig. 2 (a, b)) and the power corrections from the operators with d = 3, 
< iptp >, d = 4, m s < ijjij) >, d = 5, < ipip > (fig. 2(c, d, e)) and the 
photon interaction with a soft quark line (fig. 2f). In calculating the bare- 
loop and nonperturbative correction contributions, we first write the Lorentz 
decomposition of the correlator as: 

H-itvip, q) = ie^ aa e a q a Ui + \q tl e v - e,j,q v )Il 2 (20) 

and the dispersion representation (Cutkosky method) for the coefficients of 
corresponding Lorentz structures appearing in the 11^ (p, q) as follows: 

ds — : — — h subtraction terms (21) 
s — Q 
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Figure 2: Feynmen diagrams for bare-loop (a, b), power corrections from the 
operators with 3, 4 and 5 dimensions (c, d, e) and propagation of the soft 
quark in electromagnetic field (f) 

where pi^is, p 2 ) are spectral density corresponding to two structures in H^ip, q) 
and subtraction terms stand for corrections. To calculate pi,2(s,p 2 ), we con- 
sider Feynmen diagrams in Fig. 2 (a, b). For instance, for the contribution 
of diagram (a) we get 

d A k 



Tfu, cN c Q s J 



(2tt)< 



k 2_ m 2 W ^\ Q + k )2_ m 2 ?( p+k y_ m f»\} 

(22) 

With the help of the above equations, we obtain the following expressions 

corresponding to the coefficients of the structures ie^ uct(T e a q a and [q^e v —e^q v \. 

—eN c Q s r 1 r l 
Ti = — c s { / dxx / dy[m c m s - m 2 s xy + p 2 (X 1 y + X 3 ) - (p.q)(X 2 y 2 
An z Jo Jo 

— Ax 2 xy + 2x 2 xy 2 — Axxy)} \ dae~ a }, 

Jo 
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T 2 = ^tt%{ / dxx f dy[m c m s - m 2 s xy + p 2 (X 1 y - X 3 ) - (p.q)X 3 (2y 

lbTC z JO JO 

+ Axy)] / dae~ aA } (23) 
Jo 

where x(y) = l—x(y), A = —p 2 xxy — Q 2 xxy+m 2 x+m 2 x, Xi = x — 2x 2 +x 3 , 
X 2 = x 2 - 2x 3 , X 3 = -x + x 2 

In this calculation, we have also used the exponential representation for 
the denominator as: 

— = — / daa n ~ l e- a * (24) 

A™ (n - 1)! Jo 

Next, we apply the double Borel operator B(Mf)B(M 2 ) on T(i j2 ) and we get: 

4"7T 2 Oi + 0~ 2 JO X (Ji + <7 2 

+ p 2 {{x 3 - 2x 2 )^^ — + x 2 - x^ 1 — ) - {p.q){x 2 - 2x 3 + {3x 2 

(Ti + (72 01 + 0"2 

2 

- 4a; 3 )- — ^— - + (8x 3 -2x 2 -Ax)^ 1 — )], 

eN c Q s o x o 2 /•S.. 1 .ir m ?^v. 1+ ^r 2 o 2 



^ = eN c Q s Q 1 Q 2 /"I dx ^ e ^(mlx+m^)(a l+ a 2 )l 

16n 2 (Ji + a 2 Jo x 



167T 2 (Ti + a 2 Jo x <Ji + a 2 

+ p 2 {{x + x 3 -2x 2 ) — — x 2 + x) - (p.q)(-2x + 2x 2 " ° 2 



0"1 + U 2 (T\ + <7 2 

2 2 

2x 3 - — + (8x 2 - Ax 3 - Ax)^± — ) + {2x 3 - 2x 2 ^ ° 1 



v cri + <7 2 ) 2 ' 'gx + g 2 ' y (o-i + cr 2 ) 2 

(25) 

where <i\ = and o 2 = j^. 

In deriving Eq. (1251) . we use the definition 

B{M 2 )e~ ap2 = 5{1 - aM 2 ) (26) 

For the determination of the spectral density, we apply the Borel transfor- 
mations to Ti and T 2 [19] and we obtain: 

Qi, 2 {s,t,p 2 ) = -B(l, ai )B(~,a 2 )^- (27) 

St S t 0\0 2 



In this step, we use the following relations: 

da 

Pi,2{s,p ) = J dt — j— 2 — , (29) 

and 

B(-, (705(7, a 2 )e-^ + ^ = 5(1 - -)S(1 - j) (30) 
s t s t 

Then, we get the following expressions for the two spectral densities, as 
follows: 

Pi(s,P 2 ) = eN A C< i S 7 — org / dx={(m c m s -m 2 s x)(s -p 2 ) 2 +p 2 [(x 3 - x 2 ) 

47T Z (S — p z ) i Jxo X 

(s - p 2 ) 2 - x{p 2 + ^ + ^)( S - p 2 )] - i(m 2 - P 2 ) [(x 2 - 2x 3 ) 

Ob Ob Zi 

(s - p 2 ) 2 - 2(3x 2 - 4x 3 )(p 2 + ^ + ^) 2 + (8x 3 - 2x 2 - 4x)(p 2 

a; x 

P2(s,p 2 ) = 6 c s — — [ dx={[m c m s - m 2 s x + p 2 (2x - 3x 2 + x 3 )] 

16n z (s — p z yJx x 

(s - p 2 ) 2 -\{m 2 s - p 2 ) [(-2a: + 2x 2 - 4x 3 )(s - p 2 ) 2 + (8x 2 - Ax 3 
-Ax)(p 2 + ™k + ^)(s -p 2 ) - 2{2x 3 - 2x 2 )(p 2 + ^ + ^) 2 ]} 

XX XX 

(31) 

where the integration region is determined by the following inequality: 

sxx — (m 2 x + m 2 x) > (32) 

Similar to above calculations for diagram (a), we have repeated the entire 
calculations for diagram (b). Finally, we get the following results for the 
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spectral densities: 

Pi(s,p 2 ) = 4 8 ^2 (g ^ p 2 )3 iQ^(-^ + 5a - 5/3)p 6 + (28 + a(-103 

+a(17 + 4a)) + 52/? - a(61 + 8a)/3 + 4(10 + a)(3 2 )p 4 s - 3(-2 
+ 16a 3 + - 31/3 2 - a 2 (75 + 32/3) + a(39 + 2/3(41 + 8/3)) + 7A 2 
)p 2 s 2 + a(l - 86a 2 + /3 - 86/3 2 + a(67 + 172/3) + 18A 2 )s 3 } + 6{ 
p 6 + 2m c m s (s - p 2 ) 2 + (2 - 11a + 4/3)p 4 s(-l + 26a 2 + 2/3 2 

-2a(7 + ll/3))pV - a(l + 18a 2 + 2/3(1 + (3) - 4a(4 + 5/3))s 3 } 

1 + A + g c[A | 3( _ 3 + 5f3 _ 5a)p6 + (2g + p(_ m 

1 + a — p + A 

+/3(17 + 4/3)) + 52a - /3(61 + 8/3)a + 4(10 + /3)a 2 )p 4 s - 3(-2 

+ 16/3 3 + a - 31a 2 - /3 2 (75 + 32a) + /3(39 + 2a(41 + 8a)) + 7A 2 

)p 2 s 2 + 13(1 - 86/3 2 + a - 86a 2 + /3(67 + 172a) + 18A 2 )s 3 } + 6{ 

p % + 2m c m s (s - p 2 ) 2 + (2 - 11/3 + 4a)p 4 s(-l + 26/3 2 + 2a 2 

-2/3(7 + lla))pV - /3(1 + 18/3 2 + 2a(l + a) - 4/3(4 + 5a))s 3 } 
. 1 + P — a — A... 
N+^-a + A ]} ' 

P^P 2 ) = la C 2< 1 2 , 3 {Qs[A{-(29 + 2a 2 + /3(17 + 2/3)-a(19 + 4/3))p 6 
967T Z (s — p A y 

+ (4a 3 + a(5 + /3)(5 + 4/3) - a 2 (23 + 8/3) + 6(8 + /3 - /3 2 ))p 4 s 

+ (-23 + 2a 3 - 5/3(1 + 4/3) + a(-9 + 2/3(21 + /3)) + a 2 (-30 
48 

+ T2 — TTjL + a ( 13 + a (- 5 + 22a ) + ^ 

A^ — (1 + a — p) A 

-44a/3 + 22/3 2 )s 3 } + 6{-2p 6 + m c m s (s - p 2 ) 2 + (4 + 5a)p 4 s 

-2(1 + a + 4a 2 - 2a/3)pV + a(l + 4a(a - /3))s 3 } 

ln ]~ a + ( l~ X A + g c [A{-(29 + 2/3 2 + a(17 + 2a) - /3(19 + 4a)) 
1 — a + /3 + A 

p 6 + (4/3 3 + /3(5 + a) (5 + 4a) - /3 2 (23 + 8a) + 6(8 + a - a 2 ))p A s 
+ (-23 + 2/3 3 - 5a(l + 4a) + /3(-9 + 2a(21 + a)) + /3 2 (-30 
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48 

~ Aa + A2 _ (1 + / g_ a)2 ))P V + /5( 13 + ^(-5 + 22 /?) + « 
-44/3a + 22a 2 )s 3 } + 6{-2p 6 + m c m s (s - p 2 ) 2 + (4 + 5/3)j9 4 s 

-2(1 + (3 + 4(3 2 - 2(3a)p 2 s 2 + (3(1 + A(3((3 - a))s 3 } 
1-/3 + a -A 

o - - - J } (33) 

1 — (3 + a + A 



where a = ^ and A = VI + « 2 + /3 2 - 2a - 2(3 - 2a(3. 

The next step is to calculate contributions coming from the power cor- 
rection terms. After standard but lengthy calculations for the contributions 
of the diagrams (c, d, e ), we get: 

tt / \(cde) m c m s ,fn 2 1 m 2 m 2 r 2m 3 

ni(p,?) M ' e) = W .<ss>-—<ss>[^ + — + -^} + ^<ss>[-^ 

2mf m c m c 2m 3 m\ ,2ml 2m 3 

+ ^r4 + + + -r4 - -r < ss > r " 



ry ^2 ,yi 2 qyf^qy qyf qy'ii ryt ryf^t ryl^try tyf^ty^i 

2m c m c 2m 3 
H H H , 

2^^2^i ryf ry2i ryf qy'i 

m c m H r m 2 m 2 m 2 r 2m 3 



c 



ry! ry ^ ryf2iry rylqyi ^ ryftjry 

2ml m c fn c 2m 3 m 2 , r 2m 3 2m 3 

+ ^r4 + + -t4 + —4 + -7- < ss > r ' 



ryf t 2iiy'2i ryt^ry eyt ' qy*2i ry! ^"3 ryt^ry ryf^ry^ 

2m c m c 2m 3 

+ ^2 + ;vf ] (34) 

where r 2 = p 2 — m 2 and r' 2 = Q 2 — m 2 c . Finally, we calculate the contribution 
of diagram (f). For the calculation of this diagram corresponding to the 
propagation of the soft quark in the external electromagnetic field, we use 
the light-cone expansion for the non-local operators. After contracting the 
c quark lines in 

IWp, q)=if d A xe^ x < 7 (g) | T{s(0) 7 ^(0)c(x) 7l/ (1 - l 5 )s(x)} | > 

(35) 
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we obtain 

d 4 X ( 27r )2 m 2_ k 2 < I + m c)7^(l - 75)s} I > 



2 < 7(.9j I «T/x(.^ + m c )j v {l - 75JS) I U 

(36) 

To calculate the matrix element appearing in the above equation, we use the 
following identities: 



i 

llfloTIv = 9mlu + gaulu - g^Ja + itfMvaalalb (37) 



and photon distribution amplitudes (DA's) for twist 2, 3 and 4 

< I sivS I > = -%/37 / dui) {V \u)x e F ev (ux) 

I Jo 

< 7(9) I S7 a 7 5 S I > = j^/ 37 / du^ A) {u)x 6 F ea (ux) 

4 jo 

< 7(g) I sa a/ 3S I > = Q s <ss > du(f){u)F a p{ux) 

J 

/" dux 1 A{u)F { a a{ux) 

16 jo 



Q s < ss > 
8 



duB{u)x p (xpF ap {ux) — x a Fj3 P (ux)) 

(38) 



where F^ v is the field strength tensor of the electromagnetic field, which is 
defined by 

F^(x) = i(e v qn - e (t q u )e tqx 

F^(x) = -e^apFa^x) (39) 

The asymptotic expression for the photon wave function <p{u) in terms of 
magnetic susceptibility of the quark condensate, x(/-0> a ^ a re-normalization 
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scale (/i = 1 GeV 2 



) is defined by: 



4>(u) = x(v)u(l - u) 



(40) 



Other functions used in Eq. (1381) are defined by [201 EI] 

^ {V \u) = -20w(l -u)(2u- 1) + —(^-3^)^(1 -u)(2u- 1) 

16 ' ' 

x(7(2«- l) 2 - 3), 
^(u) = (l-(2 M -l) 2 )(5(2 M -l) 2 -l)x^(l + ^-A^ ); 

A{u) = 40u(l-u)(3k-k + + l) + 8(& -3&)[u(l-u)(2 + l3u(l-u)) 
+2m 3 (10 - 15m + 6u 2 )lnu + 2(1 - w) 3 (10 - 15(1 - u) + 6(1 
-u 2 ))ln(l -«)], 

= 40 rrfa(4-a)(l + 3A; + )[— + -(2a-l) 2 ] (41) 



where fc, £2, £2" an d /s7 are constants (see [201 EI])- Using the above 
relations in Eq. fl36l) . we obtain: 




n M ^(p,g) 



du%l) {V) (u)x a F a ^ - -^-f^g^k a j duip {V) (u)x a F aa + -j-fa-y 
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After performing integration over x and k, the following results corresponding 
to the coefficients of two invariant structures ie lluaa e a q <T and [q^Sy — e^] are 
obtained as follows: 

m c Q s <ss> f l 1 1 f 1 10 



. m m c y s < ss > r r L , i i , , w 

ni(p,g) u; = - / du^){u)— — / duA{u){ 

2 Jo p 2 — mj 16 jo (p 2 — m 2 ) 



e 2 ) 2 



3m? 



-I )1 

(p 2 — m 2 Y ' 



n 2 (p,g) {/) = -^/ 37 ^V (A) N(^-^ + , 2 2ml 2 , 2 ) - m c Q s <ss> 
A Jo p 2 — (jr — m 2 ) 2 

J 1 , , 1 1 f 1 , w 10 8m^ 

/ du<j){u)— 2~T^ duA[u){— ^ + r~2 2^3 ) 

jo p — m 2 lo jo (p — m 2 ) 2 (p 2 — m^Y 

(43) 

These results are the final results of the QCD part (OPE expression) of the 
correlator. The next step is to equate Eq. fl20l and Eq. f|T9|) (the physical 
or phenomenological side of the correlation function) and perform the Borel 
transformation, with respect to the momentum of D* meson (p 2 — > M|), 
in order to suppress the contributions of higher states and continuum. We 
obtain the following sum rules for the transition form factors, namely: 

f va\Q 2 ) = ^e^tB{ f ds^^ + II^ +e+/ } (44) 

JD* J(m c +m s ) 2 S — Q 

where V and A are correspond to 1 and 2 in r. h. s., respectively. In Eq. (144")) . 
in order to subtract the contributions of the higher states and the continuum, 
quark-hadron duality assumption is used, i.e. it is assumed that 

p higher states^ = p OPE^Q^ _ ^ ^ 

In the calculations, the following rule for the Borel transformation is used: 

1 n e^B 

B {p 2 -sY = ^ r(n)(M|)--i ( 46 ) 
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4 QCD sum rules for the form factors in- 
duced by electromagnetic penguin 

The effective Hamiltonian for the b — > 57 transition can be written as follows: 

H = -£^V tb V; s C 7 ^)sa^m b ^^ + mJ-^]bF^ (47) 

In order to obtain the transition amplitude, we need to calculate the following 
matrix element: 

<D* s \sa (lv {l± l5 )q v b\B c > (48) 

At q 2 = 0, we can write this matrix element in terms of the two gauge 
invariant form factors Ti(0) and T^O) 

<D* s (p,e^)\sa^b\B c (Q)> = ie mPX e^ VQ A Ti(0), 

< DfoeW) I sa^lsb \ B C {Q) > = [(m| a - m%)ef^ 

- (e^\q)(p + Q),)T 2 (0) (49) 

Using the relation 

<V75 = —^ £ ^P a ( 50 ) 

one can immediately obtain that T 2 (0) = |Ti(0). Then, we need to calculate 
only the form factor Ti(0). For this aim, we define the following three point 
correlation function: 

n„ a = -J d*xd*ye« Qx -^ < I T{c(y)j a s(yX0)a^q%0)c(x)i l5 b(x)} 

(51) 

where crf a s and ci^^b are the interpolating currents of D* and B c mesons, 
respectively. 

After inserting the hadrons full set with quantum numbers of correspond- 
ing interpolating currents (see also [12]), we obtain the following expression 
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for the phenomenological part of the correlation function: 

= % ^ Bc7nBc — fe^^ 2^^ q/ 3aP /3 <5 A Ti(0) + excited states. 

(m b + m c ) (p 2 - m 2 D ,)(Q 2 - m 2 B J ^ ^ 



(52) 



For the calculation of the QCD part, we write the Lorentz structure in the 
above correlator as: 

n MQ = ^ A /Q A n(p 2 ,Q 2 ) (53) 

where 

Tl per (p 2 ,Q 2 ) = --^ fdsds'— subtraction terms (54) 

[ZTl) 2 J [S — Q 2 ){S' — p 2 ) 

The standard calculations lead to the following result for the pertubative 
part (bare- loop diagram): 

ff er (s', s) = 4N c [m b m c (Ai + M + J ) - m 2 ^ - 2A 3 ] (55) 

where 

M = jg^[AS + ml-rti)-\(5 + W + m% 



m 2 \m\s' + (m 2 — s'){s — s')] 
oTTZ a >\2 ' 



A 2 = jz —[-(S + s')(m 2 b -m 2 c -s) + S(s' + m 2 c )}, 

(s — s ■ Y I 

A 3 



2(§-s'y 

= -ijibj (56) 

The integration regions over s and s' are obtained from the following inequal- 
ities: 

m 2 c <s'< Sq, 



s 



s'm 2 



m 2 . - v 



7 < s < s~ (57) 
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Figure 3: Feynmen diagrams for gluon corrections 

The quark condensate terms give zero contribution after applying the double 
Borel transformation, with respect to the p 2 (p 2 — > Mf ) and Q 2 (Q 2 — > 
Mf). Only the gluon condensates can contribute to the form factor. Fig. 
3 shows such type of diagrams. After lengthy calculations for the gluon 
condensates contribution and equating the phenomenological and QCD parts 
and applying double Borel transformation with respect to the p 2 and Q 2 , we 
find the following expression for the form factor Ti(0): 

71(0) = ~ , A 2 , + , ^ [ / d~sds'ff-{s\ i)e'M 
(2tt) 2 fBjn 2 B J D1 m D * J 

+ M 2 M 2 < —G 2 > C G 2] (58) 

7T 

where C G 2 is the Wilson coefficient of the gluon condensate and we thus have 
(see Fig. 3): 

Cq2 = Cq2 + Cq2 + Cq2 + Cq2 + Cq2 + Cq 2 (^9) 
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The explicit expressions for C G2 are given below as follows: 

C£ 2 = 96m 6 {K(7 [l,3,l]+m^o[l,4,l]+/i[l,3,l]+m 6 2 7 1 [l,4,l] 
+ I 2 [l, 3, 1] - m 2 / 2 [l, 4, 1])] + m b (h[l, 3, 1] + m 2 b h[l, 4, 1] 
+ 2/ 3 [l,4,l])}, (60) 

C b G2 = 16{2/o[l,l,2]+2m 6 m c /o[l,l,3] + 2/ [2,l,l] + 3m 6 m c /o[2,l,2] 

+ 4m 2 / [2, 1, 2] + 4m b m%[2, 1, 3] + 2m 6 m c / [3, 1, 1] + 2m 2 / [3, 1, 1] 

+ 6m 6 m^o[3,l,2] + 2m^o[3,l,2] + 2m 6 m^ [3,l,3] + 27i[l,l,2] 

- 2m 2 /! [1, 1,3] + 2m b m c h [1, 1, 3] + 2/i [2, 1,1]- mg/i [2, 1, 2] 

+ m b m c h[2, 1,2] + 4m 2 /i[2, 1,2] -4m^/ 1 [2, 1,3] + 4m 6 m^i[2, 1, 3] 

- 2mgii [3, 1,1] + 2m^i [3, 1,1] - 6m 2 b m 2 c h [3, 1,2] + Am b m 3 c h [3, 1,2] 

+ 2m 4 / 1 [3,l,2] -2m^/ 1 [3,l,3] + 2m 6 m^/ 1 [3, 1, 3] + 2m 6 m c / 2 [l, 1, 3] 

+ I 2 [2, 1, 1] + m b m c I 2 [2, 1, 2] + 4m b m3/ 2 [2, 1, 3] - 4m 2 / 2 [3, 1, 1] 

+ 4m b m%[3, 1, 2] + 2m 6 m c 5 / 2 [3, 1, 3] - 4/ 3 [l, 1, 3] - 4/ 3 [2, 1, 2] 

- 8m%[2, 1, 3] - 8/ 3 [3, 1, 1] - 16m 2 / 3 [3, 1, 2] - 4m 4 / 3 [3, 1, 3]} 

32M2 dJl 2 {M2 [/o[2, 2] + 2mbmJo & lj 3 1 + 2m ^ J o[ 3 ' !> 2 ] 
+ m c 2 / [3, 1, 2] + 2m b m% [3, 1, 3] + h [2, 1,2]- 2m 2 /! [2, 1, 3] 

+ 2m b m c h[2, 1, 3] - 2m 2 /! [3, 1, 2] + m b m c h[3, 1, 2] + m 2 c h[3, 1, 2] 

- 2m 2 m 2 / 1 [3, 1, 3] + 2m b m 3 c h[3, 1, 3] - J 2 [2, 1, 2] + 2m b m c / 2 [2, 1, 3] 

- 2/ 2 [3, 1, 1] + m b m c I 2 [3, 1, 2] - m 2 / 2 [3, 1, 2] + 2m b m 3 c I 2 [3, 1, 3] 

- 4J 3 [2, 1, 3] - 6/ 3 [3, 1, 2] - 4m 2 / 3 [3, 1,3]]} 

- 32M 2 4 (J^) 2 {M 2 4 K/ 1 [3, 1, 3] + 2J 2 [3, 1, 2] - m b m c (I [3, 1, 3] 

+ J 1 [3,l,3] + / 2 [3,l,3]) + 2/ 3 [3,l,3]]} , (61) 

C C Q2 = 96m c {-(m 2 m c / 1 [4,l,l]) + m 6 (/ [3,l,l]+m 2 / [4,l,l] 
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+ 7^3, 1, 1] + m 2 c h[A, 1, 1] + J 2 [3, 1, 1] + m 2 7 2 [4, 1, 1]) - 2m c 7 3 [4, 1, 1]}, 

(62) 



C^ 2 = -32m 6 m c {/o[2,l,2] + /o[3,l,l]+m2/ [3,l,2] + /i[2,l,2] 
+ 7^3, 1, 2] + 7 2 [2, 1, 2] + m*7 2 [3, 1,2]- 47 3 [3, 1, 2] 
- M 2 ^[M 2 (7 [3, 1, 2] + 7^3, 1,2] + J 2 [3, 1, 2])]} + lQ{m 2 b h% 1, 2] 
+ J 2 [2, 1, 1] - m b m c (/ [2, 1, 2] + 7^2, 1, 2] + J 2 [2, 1, 2]) + 4/ 3 [2, 1, 2]}, 

(63) 

C e G2 = 16{27 [l,l,2] + 7 [l,2,l] + 2m 2 7 [l,2,2]+m fe m c 7 [l,2,2] 

+ m 2 7 [2, 1, 2] + m b m c I [2, 2, 1] + m 2 / [2, 2, 1] + m 2 m 2 / [2, 2, 2] 
+ m b m%% 2, 2] + h[l, 1, 2] + 7^1, 2, 1] + mg/i[l, 2, 2] 
+ m 6 7n c /i[l,2,2] - A [2, 1,1] - 2m 6 m c 7 1 [2, 1, 2] -mg/^2,2,1] 
+ m 2 ^, 2, 1] - 2mgm c J 1 [2, 2, 2] + m b mlh% 2, 2] + 7 2 [1, 1, 2] 
+ +7 2 [1, 2, 1] + m 2 7 2 [l, 2, 2] + m 6 m c 7 2 [l, 2, 2] - 2m fe m c 7 2 [2, 2, 1] 
+ m 2 7 2 [2, 2, 1] + m 6 m;?7 2 [2, 2, 2] - 27 3 [2, 1, 2] - 47 3 [2, 2, 1] 

- 2m 2 7 3 [2, 2, 2] - 4m 6 m c 7 3 [2, 2, 2]} 

" 16M i^2{ M i[ / 2[2,2,l] + 27 3 [2,2,2]]} 

+ 16M 2 ^{M 2 2 [-m b m c 7 [2, 2, 2] + 7 X [2, 1, 2] + m 2 7x[2, 2, 2] 

- m^mJip, 2, 2] + 7 2 [2, 1,2] + 7 2 [2, 2, 1] + m 2 7 2 [2, 2, 2] 

- m fe m c 7 2 [2, 2, 2] + 27 3 [2, 2, 2]]} + \Q{h[l, 1, 2] + m^fl, 2, 2] 
+ 7 2 [1, 2, 1] - m b m c (I [l, 2, 2] + 7^1, 2, 2] + 7 2 [1, 2, 2]) 

+ 27 3 [1,2,2]}, (64) 

C f G2 = 16{27 [l,2,l] + 27o[2,l,l] + 2m 2 7 [2,2,l] 
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+ Qm b m c I [2, 2, 1] + 2m 2 / [2, 2, 1] + 2h[l, 2, 1] 

5/i [2,1,1]- 5m=/i [2, 2,1] + 6m 6 m c /i [2, 2, 1] + 2™^ [2, 2, 
+ 2/ 2 [l, 2, 1] - J 2 [2, 1,1]- m 2 / 2 [2, 2, 1] + 6m b m c I 2 [2, 2, 1] 
+ 2m 2 / 2 [2, 2, 1] - 14/ 3 [2, 2, 1]} - 32M 2 -^{M 2 [/ [2, 2, 1] 



+ 1^2,2,1] + I 2 [2,2,l]]} 



dMl 



and the for explicit form of the Ii[a, b, c], we obtain: 

(_1 y+b+c 



I [a,b,c] 



h[a,b,c] 



I 2 [a,b,c] 



I 3 [a,b,c] = 



— (Mh 2 - a ~ b (Mt) 2 - a - c 

i6tt 2 r(a)r(6)r(c) v 1J v 21 

xU (a + b + c- A, 1- c-b) , 

/_i \a+6+c+l 

i6^r(V(^r(c) (M?)2 "'" t(M|)3 "°" C 

xW (a + 6 + C-5, 1 — c — 6) , 

— \ ; - — — (M 2 ) 3 - a ~ b (M 2 ) 2 - a - c 
i67r 2 r(a)r(6)r(c) v i; v 2; 

xU (a + b + c — 5,1 — c — b) , 

— \ ' — — (M 2 f- a -\Mlf- a - c 

32tt 2 r(a)r(6)r( c ) v i; V 2; 
xW (a + 6 + c-6,2-c-&) 



The function Uo(i,j), also, is given by: 

Uo(i,j) = / dy(y + M 2 + M 2 )V exp 
Jo 

where 



B - B lV 



Bo 



B 1 = 



ml 



M 2 



1 



M\Ml 

ml 
Ml Ml 



Ml + Mf 

M 2 m 2 + M|(m 2 + m 2 ) 



21 



5 Numerical analysis 



In this section, we present our numerical analysis for the form factors. From 
the sum rule expressions of these form factors, we see that the condensates, 
leptonic decay constants of B c and D* mesons, continuum thresholds s , s 
and s' , the relevant parameters in photon distribution amplitudes (DA's) 
and Borel parameters Mj,, M\ and Mf are the main input parameters. In 
further numerical analysis, we choose the value of the condensates at a fixed 
renormalization scale of about 1 GeV. The values of the condensates are[22j: 
< | M= i G ev>= -(240 ± 10 MeV) 3 , < ss >= (0.8 ± 0.2) < ^0 > and 
ttIq = 0.8 GeV 2 . The quark and mesons masses are taken to be m c (/i = 
m c ) = 1.275 ± 0.015 GeV, m s (l GeV) ~ 142 MeV [23] , m b = (4.7 ± 
0.1) GeV [22] , m D * = 2.112 GeV and m Bc = 6.258 GeV. For the values 
of the leptonic decay constants of Be and D* s mesons, we use the results 
obtained from the two-point QCD analysis: fg c = 0.35GeV [2SJI2ZII2E] an d 
fr>* = 266 ± 32 MeV^ZQ. The relevant parameters in photon distribution 
amplitudes (DA's) are taken to be x — 3.15 ± 0.3GeV~ 2 , n = 0.2, k + = 
0, Ci = 0.4, Ci + = 0, C2 = 0.3, C 2 + = 0, / 37 = -(4 ± 2) x 10" 3 GW 2 , w* = 
-2.1 ± 1.0, = 3.8 ± 1.8 [201 IZH [25]. The threshold parameters are also 
determined from the two-point QCD sum rules: So = 8 GeV 2 , sq = 45 GeV 2 , 
s' = 8 GeV 2 [21 [MIO- The Borel parameters M 2 B , M 2 and M| are auxiliary 
quantities and, therefore the results of physical quantities should not depend 
on them. In the QCD sum rule method, OPE is truncated at finite order, 
leaving a residual dependence on the Borel parameters. For this reason, the 
working regions for the Borel parameters should be chosen such that in these 
regions the form factors are practically independent of them. The working 
regions for the Borel parameters Mj,, Mf and M| can be determined on 
the condition that, on the one side, the continuum contribution should be 
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small, and on the other side, the contribution of the operator with the highest 
dimension should be small. As a result of the above-mentioned requirements, 
the working regions for this transition are obtained to be: 

4 GeV 2 < M| < 10 GeV 2 , lOGeV 2 < M 2 < 25GeV 2 , 

4GeV 2 < M 2 < lOGeV 2 . (69) 

Now, by calculating the total decay widths and taking | V cs |= 0.957 ± 
0.017 , | V cb |= 0.0416 ± 0.0006, | V tb \= 0.77±gi|, | V ts \= (40.6 ± 2.7) x 
10~ 3 [30J, < ^G 2 >= 0.012GW 4 [29], C 7 (/i = m c ) = -0.0068 - 0.02z 
[13] and tb c = 0.52 x 10~ 12 s [31j, we obtain the numerical results of the 
electromagnetic penguin(EP), weak annihilation (WA) and total branching 
ratios for this decay as follows: 



B {EP) {B C -> D* sl ) = 3.468 x 10~ 6 
B^ WA \B C -> D* a i) = 1.557 x 10~ 5 
B {Total) (B c -> D* 7 ) = 2.462 x 10" 5 (70) 

From the above results, we see that the weak annihilation contribution to the 
total branching ratio is about 4.48 times greater than that of the electromag- 
netic penguin diagram. Here, it is observed that the difference between the 
total branching ratio with sum of the weak annihilation and electromagnetic 
penguin branching ratios comes from the cross term in total decay width. 
Also our result for the total branching ratio shows that the B c — > D*^ decay 
can be measured at LHC 

Now, we compare our results of the B c — > D*^ to the results of the per- 
turbative QCD [TH], relativistic independent quark model [IB], pertubative 
QCD in standard model (SM (PQCD)) |T7] , multi scale walking techni- 
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color (MWTCM) [TTj and topcolor assisted MWTCM (TAMWTCM) [T7] 
for tb c = 0.52 x 10~ 12 s as shown in Table (1). 





B^(B C - D* sl ) 


B WA (B C D*j) 


B T ° tal (B c -> L> s *7) 


Present study 


3.468 x 10-° 


1.557 x 10" 5 


2.462 x 10- 5 


PQCD 


3.70 x 10-° 


4.94 x 10-° 


1.14 x 10~ 5 


RIQM 


2.40 x 10~ 5 


4.51 x 10~ b 


1.39 x 10" 4 


MWTCM 


(0.68- 3.42)10- 4 


(0.74- 0.81) lO^ 


(0.74- 3.57) x 10~ 4 


TAMWTCM 


(5.18- 7.23)10" Y 


(7.24- 8.13)10" Y 


(1.78- 9.95) x 10~ d 


SM(PQCD) 


1.73 x 10~ y 


5.89 x 10~ y 


7.83 x 10~ y 



Table 1: Comparison of the branching ratio for B c — > D*^y decay based on 
the r Bc = 0.52 x 10~ 12 s. 

Looking at this table, it is seen that there is a good agreement between 
the present study and the PQCD [15], in order of magnitude for the total 
branching ratio. However, our result is approximately one order of magnitude 
less than that of the RIQM and MWTCM. Also, it is one order of magni- 
tude and two orders of magnitude greater than that of the TAMWTCM and 
SM(PQCD) [17J, respectively. The ratio of B WA /B EP for the present work, 
PQCD [IS], RIQM, SM(PQCD) [TF], TAMWTCM, MWTCM are 4.48, 1.34, 
1.9, 3.4, 1.23 and 0.01, respectively. As a result of the above discussions, 
we can say that in the QCD sum rules (present study), relativistic inde- 
pendent quark model, perturbative QCD and TAMWTCM approaches, the 
weak annihilation contribution to the total branching ratio dominates the 
contribution coming from the electromagnetic penguin diagram, but this is 
not true only for the MWTCM approach. The presence of the pseudo Gold- 
stone bosons in the MWTCM leads to a discrepancy between this model and 
the other two models in [T7] (for more details see[T7]) and a part of inconsis- 
tency in the results of the different methods may be related to the different 
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magnitudes of the input parameters, getting from different references; e.g. 
we use m c (n = m c ) = 1.275 ± 0.015 GeV for the c quark masses while the 
authors of [17] use m c = 1.6 GeV and also to the nature of the methods and 
their accuracy. 

In this step, for the analysis of B c — > D*7, in the entire calculations we 
replace the s quark with the d quark. Making mo* ->mfl., fo* — > /d*, V ts — > 
Vtd, -> Vrf changes and taking m D * = 2.010 Gel/, / D * = 0.23 ±0.02 Gel/ 
[32], Vcd = 0.230 ± 0.011, V td = (7.4 ± 0.8) x lO" 3 [30J and m d = 5 MeV we 
obtain the numerical results as below: 

B( EP \B C -> D* 7 ) = 1.151 x 10~ 7 
B (WM) (5 C -> L>* 7 ) = 2.162 x 10~ 6 
B {Total) (B c -> L>* 7 ) = 2.786 x 10~ 6 (71) 

These results also enhance the importance of the weak annihilation contri- 
bution to the total branching ratio in comparing with the electromagnetic 
penguin diagram ones for the B c — > Z}* 7 . Finally, we compare our results to 
the relativistic independent quark model (RIQM) [16J for tb c = 0.52 x 10~ 12 s 
in Table (2). 

From the Table 2, it is also seen a good agreement in the order of magnitude 





B^(B C - D* 7 ) 


B WA (B C -> D* 7 ) 


B rotal (B c -> L>* 7 ) 


Present study 


1.151 x 10" Y 


2.162 x 10 _e 


2.786 x 10" B 


RIQM [16J 


5.70 x 10- y 


1.33 x lO" 5 


3.64 x 10" 6 



Table 2: Comparison of the branching ratio for B c — > D*^y decay based on 
the r Bc = 0.52 x 10~ 12 s. 

between the present study and the relativistic independent quark model. 
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In conclusion, the present study concentrated on the radiative B c — > D*<y 
and B c — > D*^ decays in the framework of QCD sum rules. The form factors 
responsible for these decays were calculated. The branching ratio for this 
decays were estimated. The results show that the B c — > _D*7 case can be 
measured at LHC in the near future. 
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